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I. INTRODUCTION 



In the recent past new approaches to studying the low-energy dynamics of QCD based on the 
idea of AdS/CFT correspondence have emerged. The idea of this correspondence was first put 
forward in 1|. The prescription of the latter consists in using a classical multidimensional string 
theory to describe a quantum four- dimensional field theory in a strong coupling regime. Namely, 
one puts into correspondence operators in the quantum field theory and multidimensional 
classical fields on the string theory side using their transformation properties with respect 
to global symmetries of the field theory and isometries of the string theory background space- 
time. Then the AdS/CFT conjecture states that the generating functional of the quantum field 
theory equals the exponent of the action of the gravity theory in which the multidimensional 
fields are confined to classical trajectories and their boundary values are set to be equal to the 
sources of the corresponding operators on the four- dimensional side: 

Z QFT [Ji{x^)] = exp(iS Gravity ci)U i{Xti) y M= Q )=Ji{x ^ ■ (1) 

This correspondence has been formulated precisely in the case of a M = 4 Super Yang-Mills 
theory on one side and type IIB string theory in the AdS§ x S 5 . The coupling gyu and the 
number of colors iV c of the field theory are related to the radius of the sphere S 5 and AdS^, 
string tension a'~ l and coupling g s t r ing as well as to the 't Hooft constant A': 

R A 

-^2 = 4:iig string N c = 9y M N c = A'. (2) 

When we move from the M = 4 Super Yang-Mills theory towards the real-life QCD decreas- 
ing the number of supersymmetries and destroying the initial conformal symmetry we start 
facing more and more complex theories on the gravity side of the correspondence. In order 
to explore QCD from the holographic point of view two basically different approaches have 
been devised, so-called "top-down" and "bottom- up". In the former, one starts with a string 
theory setup and works all the way down to warp the ten-dimensional geometry so that it 
reflects the dynamics of QCD on its four-dimensional boundary, e.g. j^j], jjj. In the latter one 
formulates a minimal model that holographically reproduces the symmetries and dynamics of 
QCD. Such models were proposed in [J], Q|, Q|. One of their advantages is that these simple 
five-dimensional toy-models allow us to study the basic properties of QCD without the com- 
plexity of the "top-down" approach. For instance, masses, decay rates and couplings of the 



lightest mesons as well as the Gell-Mann-Oakes-Renner relationship for the pion mass have 
been studied in {4], the linear confinement - in and the chiral symmetry breaking - in [?]]. 

However simple those models may be, they need to reproduce a set of exact equations that 
describe the low-energy dynamics of QCD - the so-called low-energy theorems. These theorems 
are derived from axial Ward identities in QCD and are formulated for two-, three- and four- 
point correlators. Their complete list can be found in Qj. We will focus on the following 



theorems in the Nf = 2 case 



9], 
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y J (A 2 + m 

J d^xd'y (Q(x)Q(y)) 

f d*xd% (P 3 (x)P (y)) = -^pL = 2(m u -m d )m [ dX 



V J m 2 . J (A 2 + m 2 ) 2 

-(m u -m d ) . (5) 

Here i,j are the adjoint flavor indices of the scalar and pseudoscalar currents S(x) and 

F TV? ~ 

P(x): m n is the pion mass; G w = 2F n B = — — - is its pseudoscalar decay constant; G n = 

m 

B 2 L 7 g 2 6 

4(m<i — m u )— — — ; Q(x) = -tiF f/U ,(x)F ,J3 '(x) is the topological charge density; m is the 

F*fl 32 ix 

quark mass, or, in the case of different quark masses, m = -{m u + m^); V is the Euclidean 
four-volume; p(A,m) is the spectral density of the Dirac operator, and are the constants of 
the next-to-leading order effective Chiral Lagrangian (see, e.g., 

The aim of this paper is to demonstrate that these theorems are satisfied in the dual 
AdS/QCD models in the chiral limit, and therefore we will focus on the leading terms in 
the expansion in powers of quark mass m of the relevant expressions. We leave the case of a 
finite mass to future studies. 

In the next section of the paper [TT] we compute the necessary two-point correlators from 
the AdS/QCD point of view, in the section ITTT1 we discuss the implications of the AdS/QCD 
models for the Chiral Lagrangian, and in the section IIVI we present a way to calculate the 



3 



aforementioned spectral density in the AdS/QCD framework. Finally, in the section IVl we 
discuss the compatibility of the low-energy theorems with the AdS/QCD models. 



A. Action and fields 



We will consider the following five-dimensional action of the AdS/QCD models 

Sbd = f d^ge'Hr {A 2 (\DX\ 2 + |^|A| 2 + - ±(F 2 + F 2 )} . (6) 

R 2 R 2 
with a metric ds 2 = — i—dz 2 + dx^dx 11 ) = gMNdx M dx N = —rjMNdx M dx N . 
z 2 z 2 

Here g = det(gMN), $ is a dilaton whose profile depends on the choice of a particular 



model. In the so-called hard-wall model 



4j = and the bulk five- dimensional space has a 



boundary at z = z m where a uniform Neumann boundary condition is imposed. In the soft-wall 



models 



121 ] the bulk space stretches over the whole range ^ z < oo and the dilaton profile 



is set to be asymptotically parabolic: 

$0) ~ \z 2 , z -> oo (7) 

so that the linear Regge trajectory is reproduced. 

We introduce two gauge fields and R^ in the adjoint representation of the SUi{Nf) and 
SUii(Nf) gauge groups respectively with curvatures Fl = dL — iLAL, Fr = dR — iRAR, and a 
scalar field X af3 with mass —3/R 2 in the bifundamental representation of SU^Nf) x SUn(Nf) 
that interacts with the gauge fields: DX = dX — iLX + iXR. 

According to the AdS/CFT prescription ([T|) the fields on the AdS boundary act as sources 
of the QCD currents: 



L a tx (x,z = ti) = source of qL{x)l^t a q L (x) 
0) = source of g j R(x)7 /J t a g R (x) 
: source of q1{x)q l3 R {x). 



R°(x,z 



lim -X af3 (x, z 

z^O z 



(8) 



In the AdS action g 5 is the 5D coupling constant defined by means of a comparison of the 

127T 2 

vector two-point correlator with the QCD sum rules 



A is the normalization 



factor for the scalar field which is fixed by comparing the pseudoscalar two-point correlator in 
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the large momentum regime with the sum rules 

LE 1 

is necessary in the soft-wall models 12] 



A 2 R 3 



Nc_ 
Air 2 ' 



Hard-wall models 



generally do not include the effective potential of the scalar field V(X) = — 5-|^| 4 although it 

R 2 



II. PSEUDOSCALAR AND SCALAR TWO-POINT CORRELATORS 



We shall first discuss the simpler case of pseudoscalar currents. 



A. Pseudoscalar currents 



As compared to the scalar five-dimensional field, the pseudoscalar field is less sensitive to the 
shape of the effective potential V(X). For that reason we will use a simpler hard- wall model 

n ri 

[4j, [7( while calculating the (PP) correlator. We will need the following part of the 5D action: 



S 5D = J d'x^Ahr (j«9X| 2 + A|X| 2 ) . 



(9) 



According to the AdS/QCD prescription ([H]), -X al3 (z,x 



that the pseudoscalar current corresponds to iq a (x) , y 5 q l3 (x) 



«-> q1(x)q^(x). This means 
2 fX-X 

z^O 



a i 



z,x 



The pseudoscalar two-point correlator is the second variation of the quadratic action with 
respect to the pseudoscalar sources. (See the appendix [S] for the details.) 



i{P i (q)P j (p))= -le 



— ?' p~ iS 5DlJp] 



5_ 5 

MM 6J 3 P (p) 



JSsd[Jp] 



__ 5r P (q) sjUp) 



SbD[Jl 



A R 5ijV(Q)5(p + q) + terms originating from the X «-> L a ,R a mixing 



(10) 



where T>(Q) is the integral of the on-shell Lagrangian with unity sources on the AdS bound- 
ary over dz, see (IA8j) : 



V{Q) 



1 4a 



2,2\ 



Q 2 \og(Q 2 e 



Here we have introduced the UV cutoff at z = e. 

The parameter a is related to the quark condensate (1A14I) 
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C = A 2 R 3 + a 



This allows us to express the pseudoscalar correlator in ( jlOl) through the condensate: 

f C A 2 R 3 1 
i (Pi(q)Pj{-q)) = Sij < — + — ^— Q 2 log(Q 2 e 2 ) V + terms from the X <-> Z£, mixing.(13) 

The X <-» L°, mixing is proportional to Q 2 , so that when Q 2 = 

i(Pi{0)Pi{0)) = Si~. (14) 

One can see that the expression has the pole — oc — ^ due to the pion exchange. 

m mi 

TV 

Let us also consider a particular Nf = 2, m u case. The difference from the previous 

consideration will arise in the expression of the bulk-to-boundary propagator and of the function 

T>(Q) = - \ + Q 2 log(Q 2 e 2 ) +0--M" 1 (c.f. (jllli ). where the matrices M = diag(m u , mA a = 
4 |_e 2 

diag(cr u , 0"^) ~ a ■ I. The second variation of the action due to pseudoscalar sources gives: 

• (P t (q)P 3 (-q)} = 5 tJ ^- Q 2 log(Q 2 e 2 )5(p + q) + 2 tr(f E*»') 5(p + g) 

+ terms from the X <-» L a , mixing. (15) 

The generalization of the result (TT4T) is the following: 

C CAm 

i (Pi(0)Pj(0)) = dij— 2~r(^ 3 3 + ^'°^ 3i _ z '^i<% - iSji5 2 i), (16) 

where Am = m u — m^. 

B. Scalar currents 

The model used in the previous subsection is insufficient when dealing with the scalar cur- 
rents due to their sensitivity to the effective potential V(X). Thus, in order to calculate the 



scalar two-point correlator we will consider a soft-wall AdS/QCD model [12] in which V(X) 
arises naturally. In this model the dynamical fields have to vanish at z — > oo and the solutions 
of the equations of motion are one-parametric unlike (1A5j) . This implies that in the absence 
of the quartic term in the action the solution for the scalar field would be proportional to its 
source - the quark mass - and would vanish in the chiral limit. This would not allow us to 

(3 



disentangle the spontaneous chiral symmetry breaking from the explicit one. Introduction of 
the effective potential restores the correct chiral limit and the constant k generates the mass 
splitting between axial and vector mesons, k ~ 15 being the best fit to the radial spectra of the 
axial mesons [3]. The parameter A in (j7j), responsible for the slope of the Regge trajectory, is 
determined to be A « 0.183 GeV 2 [l2]. 
We will use the action in the form 



S 5D = j d 5 x^e~*Ahr (\dX\ 2 + A|X| 2 + ^|X| 4 ) . 

In the case when X(z) is real and proportional to the unity matrix, X(z) 
nonlinear equation of motion assumes the form: 



V{z) 



(17) 



• I, the 



d 



-*(*) 



Z° I Z° 

The AdS/QCD prescription (JSj) implies that 



3V(z) + ^V\z)) = 0. 



;is) 



V(z) ~ mz + az 3 , z — > 0. 



Eq. ffT8l) may be rewritten as an equation for the dilaton [12] : 



(19) 



dMz) 



(X 



d z V{z) 



d z V(z 

from which it follows that asymptotically 



2z 



V{z) ~ 72;, 7 = 2\ -. 



(20) 



(21) 



The freedom in the behavior of the dilaton for small values of z enables us to choose any 
form of the function V(z) provided that it satisfies the conditions (I19ll2ip . For instance, 

.,2 



V(z) = z ( m + (7 — m) tanh 



az 



7 — m 



[22) 



is a suitable choice. The dilaton will have an asymptotically parabolic profile that switches 

from -Km 2 z 2 , z — ► to Xz 2 , z — >• 00 fl^ . 

The nonlinear dependence of the classical solution V(z) on the mass m can be interpreted 

as a nonlinear form of the bulk-to-boundary propagator in this model. More exactly, the latter 

d d 2 

is be defined by the whole set of functions V(z), — — V(z), — — xV(z) etc. in the case of a 

dm dm 2 

source that is uniform (i.e. its four-momentum equals zero) and scalar (i.e. proportional to the 



identity matrix in the flavor space lN f xN f )- Instead of a variation due to a scalar source we 
will be differentiating the action with respect to m. 
The 5D on shell action ([17]) equals: 



S 5D c, = I dz^N f h*# (±V0 5D V - ±kV* 



J -0^N f A 2 R 3 ±VO d 5D V = S vol + S surf , (23) 



where C 5D = d z ( + 3 — — O d 5D = g— ^ 



z° I z° z° 

5D I ' 



e -*(*) 

' K ~~ K 5 



On shell O rd V - -£V 3 = (c.f. (fH) and 



2 

= — iV/A 2 ^* / dz e -^-V\z). (24) 

The second derivative of the on-shell action (TTTT) with respect to the mass yields: 

i (S(O)S(O)) = \ dlV ^ + \ d m V d m V + \ d m V ^ d m V. (25) 

The details of the calculation can be found in the appendix [Bl eq. (1B2| - [B6|) . 
The result given in (1B7|) is the following: 

3 . „ . 3 . . N r 



i (^(0)^(0)) = 6 ij \—\N c A + —mV\KN c A 1 + 



7T 2 167T 2 47T 2 e 2 

3 2 / 2 , /27r 2 (7e 2 /k\ N c [tf a 

+ 32^'" K * - 3 + ' 0g (— V A )~2ikh A ^>' (26) 

where numerically 

A = 0.377, = 0.977, A 2 = -1.487. (27) 



III. THE QUARTIC PION LAGRANGIAN 

The 5D effective action of the AdS/QCD model can be interpreted not only as the generating 
functional for the correlators of the QCD currents but also as a low-energy action of mesons. 
Indeed, if the five- dimensional fields on the boundary are sources of the QCD currents and these 
currents in their turn are sources of the corresponding mesons, one can Kaluza-Klein decompose 
the fields in the bulk and obtain an effective action for the modes of the 4D boundary fields. 
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These modes are proportional to the wavef unctions of the mesons with corresponding quantum 
numbers. Integrating out all the dynamics along the z axis we will obtain an effective low- 
energy action for the mesons which may be rewritten as a sum of the lowest-order terms of the 
Chiral Lagrangian. Since both Lagrangians possess the same symmetries we shall assume that 
one obtained from the AdS/QCD point of view reproduces the QCD Chiral Lagrangian. 

The NLO part of the Chiral Lagrangian that we are interested in is the following (see, e.g., 



C x NLO = L{Tt 2 {dprfdPU) + L 2 Tr (d^d u U) Tr {dPrfVU) + L 3 Tr (8^ d»Ud u U ] d v U) 
+ L a Tt(D^D^U)Tt(U^x + X^U). (28) 

We will focus on the pion dynamics and will therefore need to consider the axial fields in 
the bulk. 

Let us consider the gauge sector of the hard- wall AdS/QCD action (c.f. (jSJ)): 

S 5D = J d 5 x^ti (-\{Fl + Fl^j , (29) 

where F L = dL - iL A L, F R = dR - iR A R, DX = dX - iLX + iXR,X(z) = -v(z) ■ 

^N f xN f ,v(z) = mz + az 3 . We can redefine the action in terms of the vector and axial fields 

L = V + A, R = V — A and chose a gauge V z = A z = <9 M V^ = 0. The A^ field is divided into 

transverse and longitudinal parts: A^ = Aj_^ + <9 M 0. 

Since for the axial and pseudoscalar currents d^A^ = 2mP, then —2m<j)(z = 0, x) is the 

source of the pseudoscalar current. In its turn, the pseudoscalar current being the source of 

C 

pions, (0\P a (p)\TT b (q)) = iGT T 5 ab 5(p — q) — i— 5 ab 5(p — q), would suggest that 4> a (z = 0,x) is 

Jit 

proportional to ir a (x). We will establish the exact proportionality later. 
We can Kaluza-Klein decompose the fields V and (/>: 

v fa x ) = E^WvH nz,x) = J2r in \x)f^(z). (so) 

n n 

The normalization is the following: 

J ^/W(z)/M(^) = 5 mn . (31) 

The functions V£ (x) and (j) a ^(x) are proportional to the four- dimensional wavefunctions 
of the vector and pseudoscalar mesons respectively. We shall focus on the lowest mode a< -°) (x) 
which corresponds to the pion field 7r a (x). 

9 



The functions f(z) in their turn are the solutions of the equations of motion with the 
boundary conditions /(0) = d z f(z m ) = 0. The lowest modes f<t>(z), fv( z ) correspond to the n 
and p mesons. fv(z) has the following form: 

z 2 

f v (z) = N v zh(m p z), N v 2 = -y (l 2 (m p z m ) - I 2 (m p z m )I (m p z m )) . (32) 

After we integrate out all the dynamics along the z axis and leave only the lowest mode of 
the field 4>(x) the five-dimensional action (129]) yields: 



S. 



d 4 x tr (n 2 ■ ^d^°\x)d^°\x) + ■ [d p ^°\x), ^ (o) (x)][^0 (o) (x), <9*> (0) (x)] 
g m4> 2-md^ \x)d^°\x)]. (33) 



The quantities N n , g^, g m ^ are the corresponding integrals over z and their explicit 
expressions can be found in the appendices [D]eq. ( lD2l) .lEleq. ( 1E4I) and[D]eq. ( 1D6I) respectively 

The canonical normalization of the pion field in (1331) demands that the proportionality 
between the pion and the assume the form: 

N^\x) = F^ a {x). (34) 



There are four terms in the next-to-leading order Chiral Lagrangian (1281) that are relevant 
to us: 

Vi = Tr 2 (d p U j d^U) , 

V 2 = Tr (dprfdJJ) Tr (d^d u U) , 

V 3 = Tr (d p U j d^Ud u U^d u U) , 

Vi = Tr(D p U^D>*U)Tr(U f x + X ] U) = 8^md tl 7T a (x)d tl 7r a (x) + C(tt 4 ). (35) 



Analogously to (l5j the Skyrme-like quartic term in (1331 induced by the AdS dynamics is 
the lowest power in 7r of a linear combination of Vi. 

g^^d p 7r(x),dMx)W7r(x),d^(x)] = Y,^V l + o(7r\p 4 ), 
*■ i=i 

where L x = -^4- ^ = ~— 4 = 2Li, ^3 = ^4 = ^i- (36) 

The relation between Li,L 2l L 3 is in agreement with [2], but, unlike in [jjj], we did not 
need to consider the pint interaction. 
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From (|33D and ([35]) we obtain (see (jDTjl ): 

L 4 = ^g m ^N- 2 . (37) 
One can explicitly demonstrate that Lj are regular in the chiral limit and are oc N c . 



IV. SPECTRAL DENSITY OF THE DIRAC OPERATOR 

While both the correlation functions and the Chiral Lagrangian parameters may be obtained 
from the effective 5D action more or less straightforwardly, the spectral density of the Dirac 
operator is a more complex entity because of its fermionic nature. In order to be able to 
calculate it in dual holographic theories we have to express it through the partition function. 

Let us consider a Yang-Mills theory of gluons with N c colors and Nf flavors of quarks 
\1/ in the fundamental representation of SU(N C ) with mass m and coupling g. We define the 
Dirac operator D as follows: 

D = 7"(0 M + i<7A M ). (38) 
Let {A n } be the eigenvalues of the Dirac operator: 

for VA G {A,„} 3^ : im = M. (39) 

The spectrum is Nf times degenerate due to the global SU (Nf)y symmetry. D and 75 anticom- 
mute, so that for every eigenvalue A n and eigenvector \l/ n there is eigenvalue — A n corresponding 
to 7 5 ^ n . 

We can define the Euclidean four-volume V and the corresponding spectral density of the 
Dirac operator: 



pw = Uj2 6 ( x - x ^) • ( 4 °) 

\ n I A 

We shall introduce a regulatory mass u and rewrite an infinite sum in terms of functional 

n 

trace. The density was originally expressed in terms of a resolvent in I17J. 



-i 



v \ ^— ' 7 tiv \^o^ /j, 2 + (x- x n y 

\ n I A \ n 

lim (Tr[iD - A - z^]" 1 - Tr[iD - A + i/j,}' 1 



2txV ^0 \ I a 

d 



1 (j 

— — lim — /log Det[iD - A - i/i] + log Det [iD-X + in]) (41 ) 

2'KV ^0 OjJi \ I A 
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Using a procedure called the "replica trick" (log z = -§^z n \ n _ Q ) we will be able to rewrite 



the logarithm of a determinant in terms of a partition function with the introduction of ghost 

n 

quarks \. Similar calculations can be seen in |16| . 



N f 

p(A) = -- lim j— 5ft ( log / 11 DxiDxi exp <j - / d 4 x xS D - A - ¥)Xi \ ) ( N f X flavors) 



7T1 



Zqcd 



i d I r 

— lim — K ( log / IT DxiDxi exp 
Vi^odfj, \ J fl 

d d I f n ' Nf f 
^dfl^dn^ij II D Xi d Xi exp j- J d A x Xi{iD - A - ifi)xi 

i a a r Nf n ' Nf 

—— lim — lim — 3? / DA TT DVjDVj 1 [ Dx/^X/ 

x exp ji y rf 4 x - A - in)xi + - ™)^J - ~tr F^F^ + J^tr F^F^J . (42) 

For further simplicity we will introduce the following effective actions: 

Nf 

S QCD = -i\og Z QCD = -ilog / DA JjD*jD*j 

^ j=i 

x exp ji y rf 4 x - m)^j - ±tr F^ U F^ U + J^tr F^F^J , (43) 

SQCD+ghosts = -« log Z QC D+ghosts = -i log / FA D^jD^j ] [ DxiDXi 

J J=l 1=1 

x exp jz y d 4 x x/(^F - /i - iA)x/ + ^j(*F - m)Vj - hi F^ v + J^tr F^ v J(44) 
The expression (I42D may now be simplified: 



iSQCD + ghosti 



1 (88 d d \ 



//=n=0 



(45) 



/i=n=0 



A. Spectral density in hard wall AdS/QCD 



Let us denote the physical quarks ^>j, J = l...Nf, and the ghost quarks xii I = l---n> • Nf 
(See formula 1421) . We use the AdS/QCD prescription (jSJ) to put a five-dimensional scalar field 
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X a/3 (z,x), where a,/3 = l...{n + l)Nf, into correspondence with the scalar currents: 



Z->0 



X a \z,x) 
X a V>+ N f\z,x) 



<- M>2(x)xJ(x), a = l...N f , (3 = l...nN f , 
<- a = l...nN f , (3 = l...N f , 



z^O 



TlWx^x), « = 1-nNf, (3 = l...nN f . 



(46) 
(47) 
(48) 
(49) 



2^0 



The Lagrangian (I42p suggests that the Lagrangian of the five-dimensional model possesses 
a SU(Nf) x SU(nNf) gauge symmetry. The action of the scalar field X is (jSJ): 



S 5D = J d'x^Ahr (j«9X| 2 + ^|*| 2 ) 



(50) 



The source for X is: 
2 



X a(3 (e, x) = &i&g(m_^_m, A + i\i . . . X + iu ) = m • P* + (A + z/i) ■ P x , 



N f 



nNf 



where P* = 5 Q/ g iff a,/3 = l...Nf is a projector on the quark states ty a in flavor space, P x = 
S a /3 iff a, (3 — Nf + l...{n + l)Nf is an analogous projector on the ghost states x 13 ■ 

The function X(z,x) in the bulk is expressed in terms of the source via the bulk-to-boundary 
propagator, similar to (1A5j) : 



f d 4 q 2 

X a/3 (z,x) = J J^y e l9X IC^{z,q)-X^{e,q) 



d 4 q 



* e-^/C^z, q)6(q) [mP* + (A + ^)P*J (51) 



(c.f. dATD). 

We assume that the bulk-to-boundary has the following form: 

K m (z, q) = JC(z, g)P* 7 + K{z, q)W ar (52) 

In QCD with quarks ^ with mass m and condensate S = CNf = aNfA 2 R 3 the field X in 

1 1 

the bulk equals X(z) = - (mz + az 3 ) = -^v(z). A QCD-like theory with ghosts x with mass 

1 (5 

\ + ifi should have the field X of the same form: X(z) = - ((A + 2/1)2: + az 3 ), where a = , 

2 A 2 R 6 
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C being the condensate of ghosts. In the first approximation C = C, but this equality is 
imprecise because of the mass-dependence of the condensate. This leads us to a conclusion 
that the quantity T>(Q) (analogous to T>(Q) from subsection HT1 eq. (1111) for quarks \& where 
that quantity was defined as the integral over dz of the on-shell Lagrangian with unity sources 
on the AdS boundary) equals: 



(53) 



From Lagrangian in (142]) it follows that the current XX nas to have a complex- valued source, 
while xinX has to have none. This implies that X(e,x) has to be complex and X(e,x) — 
X\e,x) = 0. The only way to do this is to temporarily introduce another field Y instead of 
X* and to treat X and Y independently, so that in general Y^ ^ X, in fact, Y — X . 

The action may be expressed in terms of the sources: 



5D 



A 2 R 3 tr 



d A q 



£>(Q)P* + V(Q)F* 



X [mP* + (A + ifi)F x ] 5(q) x [mP* + (A + ifi)¥ x ] 5(q) 

= A 2 R 3 5(0) (m 2 N f V(0) + (A + tfi) 2 nN f V(0)) = N f V ( mC + (A + ifi)nC ) . (54) 



To calculate the spectral density we will use the formula (j45p . Since S^d is a linear function 



of fj, ■ n, the first term -§^S-^S 



0. Hence, 



n=0 



H=n=0 



— U NfV 
nV 



d ~ 

-C + i(X + ifi)—C 



(55) 



Let us denote C{m) = ^^L C , nm™ for 9?(m) > 0. The expression (1551) yields: 

N °° 11°° 
p(A) = f - V C n (n + l)X n = —S - - V £ n (n + 1)A" (3fc(A) > 0), 



(56) 



n=0 



n=l 



where S(m) = ^2 S n m n for 9£(m) > 0, S = S| m=0 . 

n=0 

Making use of the formula from 13] that describes the mass-dependence of the quark con- 
densate E(m) 

E(m) = So (l- 3m '^^ + ...), (57) 

we obtain: 



P(A) = -^S (l - ^p^|A| - i] ^|A|log|A/^| ) . (5? 
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As one may notice, firstly, the result satisfies the Casher-Banks identity 181 ] . Secondly, it 



reproduces up to the factor oc Nj — 4 a well-known formula from QCD, first derived in 



^ oc 4 (59) 
p(0) F* V ' 

for the term linear in A. 

Nevertheless, eq. (1351) does not describe the dependence of the spectral density on the mass 
m and does not include the terms A 2 and higher powers of A. 



V. LOW-ENERGY THEOREMS FROM THE HOLOGRAPHIC POINT OF VIEW 

Having obtained in sections IH — IIVI the expressions for all the necessary ingredients of the 
low-energy theorems in the AdS/QCD framework we shall determine now whether (and how 
accurately) these theorems are satisfied in the dual models. 

We will be using the following result for the two-point correlator of the topological charge 



density 



2l|: 



f d*xd*y (Q(x)Q(y)) mBF^ 

= — h o(m) = mC + o(m), (60) 



V 



that was reproduced in AdS/QCD in 20]. 



Let us consider the first low-energy theorem ([3]). One can see that all its three sides have 
a pole at m — > 0, and the pole residue in the left-hand side C (1T4|) equals that of the central 
side. Indeed, 

Gl 4C 2 C 

As for the right-hand side, from fl56|) we obtain 

, 4m 2 p(A,m) 4Cm 2 f „ 1 C 
{A z + m z ) z 7r J {A z + m z ) z m 

which agrees with the residues of the left-hand and central sides. As for the finite part of the 

left-hand and central sides, one can straightforwardly check that both (Si(0)Sj(0)) (1261) and 

B 2 (L 3 - 4L 4 + 3) (J361 E7J) are oc N c . 

In the case of the second low-energy theorem (jlj) one can also observe the similar analytical 

structure of the left- and right-hand sides at the m — > pole, the residues are again in agreement 

(see (EHEOD) 
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As for the third low-energy theorem (j3J, one can see that the left-hand side defined by ( 1161) 

(160|) . coincides with the right-hand side when m 0: 

2CmAm f „ 1 Am, ,,,, , CAm . . 

^7^— ^ " T-rC™ + f l =-^ + O 1 . (63) 



VI. DISCUSSION 

As we can see, simple hard-wall or soft-wall AdS/QCD models allow us to test their compati- 
bility with the low-energy-theorems of QCD. Indeed, using the dual description of the quantum 
field theory in the strong coupling regime we have demonstrated that the theorems are sat- 
isfied in the chiral limit when the quark and pion masses tend to zero. The zero-momentum 
two-point correlators of the pseudoscalar currents precisely coincide with the corollary of the 
theorems ( JBTl l62| 1631) . and moreover the result ( fl4l) fort the correlation function (Pj(O)Pj(O)) 
independently agrees with the field theory result (|6ip . 

Holographic models also induce low-energy dynamics of four-dimensional fields on the AdS 
boundary that can be put into correspondence with meson wavefunctions with suitable quan- 
tum numbers. This allows us to calculate the parameters of the effective low-energy Chiral 
Lagrangian in those models ( |36l 1371) . The established relation between the coefficients L\, L2 
L3 ( 1361) reproduces the result jl5| in the "top-down" Sakai-Sugimoto model Q]. 

Finally, we succeeded in expressing the spectral density of the Dirac operator in terms of 
a partition function of a theory which includes QCD and ghost quark fields (H2| H5|) . This 



allowed us to calculate the density in the AdS/QCD framework, the result (1551) agreeing with 
the Casher-Banks identity [3] and reproducing the general properties of the QCD formula [jjj] 
in the linear order in A. 

Being for the most part simplified, AdS/QCD models such as Q|, 12] and Q] need further 
improvement and sharpening. One of their major shortcomings is the significant dependence 
of the results on the particular choice of the background geometry, presence of the dilaton 
and its profile, form of the effective potentials etc. For instance, the two-point correlator of 

4j, significantly differs from the ex pres sion 



scalar currents calculated in the hard-wall model 
( f26l) and has an unphysical pole structure. On the other hand, the soft-wall model [l2j] has 
not been yet directly generalized to describe quark flavors with different masses, and therefore 
in its present form does not allow us to calculate the spectral density using the formula (j45j. 
We intend to study the dependence of the AdS/QCD results for the coefficients of the NLO 
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Chiral Lagrangian on the choice of the model (hard- wall/ soft-wall) elsewhere, although at the 
moment it seems that the hard- wall model does not reproduce the coefficient Lj. 

All the aforementioned suggests that studying the AdS/QCD models by means of testing 
them with exact equations established on the field theory side of the correspondence is quite 
a promising task and will hopefully lead to the refinement of those models. Alongside with 
the advancement of the "top-down" approach that produces more rigorous results at a cost of 
significant complexity, this should reveal us the possible structure of the ultimate dual of the 
QCD. 
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APPENDIX A: THE VARIATION OF THE ON-SHELL 5D ACTION WITH RE- 
SPECT TO PSEUDOSCALAR SOURCES 

We consider the scalar sector of the 5D action ([6]): 

S 5D = f d 5 x^A 2 tx (j«9X| 2 + JjIAT) , (Al) 

see (|2j). 

If X is proportional to the unity lN f xN f matrix, the action in components is rewritten as 
follows: 



S 5D = I d 5 x N f A 2 R 3 



d^Xd^X - {d z Xf 3 v7 





Z 3 



(A2) 



The equation of motion for the bulk-to-boundary propagator )C(z, x — y) is: 



/ i \ i 3 2 

d z -^d z 1C(z,x - y) -—nK,(z,x-y) + —JC(z,x-y) = 0, \im -JC(z,x-y) = S(x-y.) (A3) 
\Z J J z^ z a z-*o z 
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The solution in the 4D Euclidean momentum space is 

K{z, Q) = Q 2 z 2 (AKx(Qz) + Bh(Qz)) . 



(A4) 



The boundary condition implies that A 



a 



2Q : 



and the Q — > limit - that B 



mQ 3 
Thus 



'1 + 0(1)), where a is proportional to the chiral condensate C, {(fq 13 ) = S a/3 = CS a/3 . 



JC(z,Q) = 9^ Kl {Qz 



az^ 



h{Qz) 



(A5) 



2 ' m Q 

The source of the pseudoscalar current Jp \x),a, (3 = l...Nf, is related to the boundary 
value of the X field on the regulatory UV brane: 

2 
e 

allowing us to express X in terms of J: 



-X a/3 {e, x) = ijf(x) + m5 a(3 , e -> 0, 



(A6) 



X a p(z,x) 



d 4 q 
d 4 q 



iqx K(z, q) U P (q) + ^W f m5(q)S 0i 



l af3i 



e- iqx K{z, -q) -iJ l P (-q) + j2N f mS(q)5' 



where t % a p are the basis Hermitian Nj x Nf matrices, % = 0...N 2 — 1. 



If we denote the integral over z as: 



dz 



z 3 z 5 



}C(z,Qy 



K{z,Q)d z K{z,Q) 



(A7) 



^ + - + Q 2 iog(Q 2 ^ 



C m 

the action may be rewritten in terms of the pseudoscalar sources: 

S 5D [J P ] = A 2 i? 3 ^ / 

11 2 J (2tt) 4 



(A? 



-StplJp], 



V(k) iJ^(k) + ^/2N~ f mS(k)5 0a -iJp(-k) + yf2N~ f m5(k)5 0b 

(A9) 



where S^[Jp] is the action for the longitudinal component of the axial field <fi which mixes with 
The pseudoscalar two-point correlator is the second variation of the quadratic action: 



i(P i (q)P J (p)) = 

5 S 
= SJp(q)SJ P (p) 



-i e 



-iS 5D [J P ] 



s 



s 



SJ P (q)SJ p (p) 



J P =0 



S 5 d[Jp] = A 2 R 3 S ij T>(Q)5(p + q) + terms originating from S^JpfAlO) 



In order to establish the precise connection between a and the condensate we will calculate 
the latter as a variation of the 5D action due to a scalar source in the same way as we have 
calculated the pseudoscalar correlator [3] . In that case the boundary condition is the following: 

2 



-X a/3 (e,x) = Jf(x)+m8 af) 



(All) 



and the 5D action is 



S 5D [J s }=A 2 R 36 f 



d 4 k 



Variation of the action yields: 



V(k) \j a s {k) + ^2N~ f m5{k)5 0a \ \j b s {-k) + y^2N~ f m5 (k)S ob \ A12) 



Thus, 



E a/3 (g) = ^{q)t a a 



ai3 



—it a p- iS 5D[ J s] 



,iS 5 o[Js] 



J S =0 



S 



JJ«(q) S5D[Js] 



Js=0 



t a ap A 2 R 3 ^-2V(k)6 at 6(k - q)^W f m5(k)5 0j 



t%^2N~ f A 2 R 3 mV(0)5(q) = 5 al3 A 2 R 3 + a) 6{q) 



J af3 



C5 af) = 5 a pA 2 R 3 { 



/ m 



V4e 2 



a 



C = A 2 R 3 ( 



/ m 



V4e 2 



a 



(A13) 



(A14) 



This allows us to express the pseudoscalar correlator through the condensate: 



C A 2 R? 

i {Pi(q)PA-q)) = <j — + — — Q 2 \og(Q 2 e 2 ) }> + terms originating from SJJ P ]. (A15) 

m 4 



APPENDIX B: THE VARIATION OF THE ON-SHELL 5D ACTION WITH RE- 
SPECT TO SCALAR SOURCES 

In this appendix we will discuss second derivative of the on-shell action ( TT7I) with respect to 
the mass: 

i (S(O)S(O)) = \ dly ^ + \ d m V 5 -^f d m V + \d m V 6 -^ d m V. (Bl) 
The first term here 



~ dly = if d 2 m v{zf—^v{z) + v(z) e —^d 2 m v(z] 



2 m 5X 2 m z 3 z 



0, (B2) 
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(7 - mf CORh 3 ( a* 



The second term in ( l25l) 



7— m 



1 ~ -. S 2 S surf N f A 2 R\ ,„ N e 

d m V ~ 



5X 2 



*(*)/) 

<9 m y(z r - - d m V(z) 



1 1 



Ar /A ^(---^). (B3) 



The third term in (1251) is the most involved. 



no/ 



4 m 5X 2 



00 p -*(z) 



d m V = -rNfA'R'n / d« — =— V\z)d m V\z) 



N f A 2 R 3 Kj 2 I dx 



00 e -*(z(a:)) 



??2 



tanhx H (1 — tanhx) 

7 



[l — tanh x + x cosh 2 2] 2 (B4) 



<7z 2 era 2 

where in the last integral we have introduced a dimensionless variable x = , e 



The functions 



7 — nC x 7 — m 



fo(x) = tanh 2 x ( 1 — tanhx + 



x 



cosh x 



/x, 



fi(x) = 2 tanhx(l — tanhx) ^1 — tanhx + 
f2(x) = (1 — tanhx) 2 I 1 — tanhx 



x 



x 



cosh x 



cosh 2 x 
2 

fx, 



/x, 



decrease rapidly when x > 1 and we will use the z — > asymptotic behavior of the dilaton: 



Km 2 2 Km 2 (7 — m) 



la 



X = XX. 



Numerically 



A ( 



dx f (x) = 0.377, Ai = / = 0.977, A 2 = / dx / 2 (x 



-1.487. 



(B5) 



As a result, 
1 „ . . S 2 S„ 



d m V d m V = ^XN f N c A + JLmyfaNfNcA! 

0X 2 8n 2 lQn 2 1 



+ ^m z KN f N c 



2ir 2 Ce 2 [k\ N c /A 3 



N c V V 2tt 2 C V « 



A . (B6) 



We finally obtain from (I251IB2IIB31 IB4IIB6]) : 

1(^(0)5^(0)) = 5 l3 

3 



3 3 _/y 

AiV c Ao + TTr^mVAfi-iVcAi + 



8tt 



16tt 2 
2 



32tt 2 



7Jl 2 KN c \ A 2 h lo£ 



47r 2 e 2 

27r 2 Ce 2 [k\ N r /A 3 



N r V A / 2tt 2 C V « 



A \y (B7) 
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APPENDIX C: THE LOWEST KALUZA KLEIN MODE OF THE FIELD (f> 

In this section we will explicitly derive the formula for the lowest Kaluza-Klein mode = 
ftj, of the field <p by solving the corresponding equations of motion perturbatively up to the first 
order in m and m 2 oc m. 

In order to find the first-order corrections to the chiral limit of the solution we will solve 
the following system of equations of motion perturbatively: 

d z Q d z ^j + fc 2 ^(7r - 0) = (CI) 

mld z <j ) + k 2 ^d z 7r = 0. (C2) 

Here k 2 = A 2 R 2 g 2 = 3, A 2 R 3 = A 2 R 3 a = C. The small parameter of the perturbative 
solution is m and m 2 oc m. The boundary conditions are: 

tt(0) = 1, 0(0) = = 0. (C3) 

We cannot impose uniform zero boundary conditions on both functions <fi and it without making 
the equations incompatible. 

The solution in the chiral limit is: 

7r(°)(z) = const = tt(0) = 1, (C4) 

= '-ir'Wlfy)" ^ ^) 

- 2r- 1 (l/3) ^— J zK 1/3 (kaz 3 ). (C5) 

We will impose the boundary conditions (1C3|) at z = and 2 = z m upon the Green function 
G{z,z'). 

The Green function is a solution of the following equation (ICli : 

<9 2 - fc z') = 5(z - z'), g(o, z') = 0, d z g(z, z') \ z=Zm , (C6) 
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t' /-2/ 3 (*m) {K 1/3 (t')I„ 2/3 (t>) + K_ 2/3 (t')/l/3(t')} 

x <j K 1/3 (t')K_ 2/3 (t m )I 1/3 (t) + 0(t - t')h /3 {t')I_ 2/3 {t m )K l/3 {t) 



+ 6{t' - t)K 1/3 (t')I- 2/3 (t m )I 1/3 (t) >, t = kaz 3 , t' = kaz'% t m = kaz 3 m . (C7) 



z z 



t' I-2/ 3 (t m ) {K 1/3 (t')I_ 2/3 (t') + K. 2/3 (t')I 1/3 (t')} 
x S^K 1/3 {t') (K_ 2/3 (t m ) + I- 2/3 {tS) h/s(t) 

+ 9(t - t')I_ 2/3 {t m ) U/ 3 {if)K 1/3 {t) - K 1/3 {t')I 1/3 {t)\ } (C8) 
= Fo(z') zl 1/3 (kaz 3 ) + 9{z - z') ^Fr(z') zl 1/3 (kaz 3 ) + F K (z') zK 1/3 (kaz 3 )j .(C9) 



The first-order correction to ir(z) is (jC2[) : 



/c 2 Jo v 2 (it) 
and the first order correction to is fIClj) : 

= / "dro ^(2,^)^(7^ (0 (o) (w) -vr (0) H) - A; 2 cr 2 W7 3 7r (1) (u;)^ . (Cll) 







We have retained the terms oc m, m 2 in the v 2 (u) function in the expression fIClOl) for the 
sake of convergence of the integral although we exceed the necessary accuracy level. 
Substituting (JCH ESQ ICTt)]) into flCTTj) and using the form (|C9} we obtain: 



<f> W (z) = -21^(1/3) f^j j^dw g(z,w) ^2mak 2 w 2 [h/s(kaw 3 ) f 2 J^ 

+ K 1/3 (kaw 3 ) + 3ka 3 m 2 w w 3 / du — - 
J Jo v 2 {u) 



x l i -"*'""'3§" ,f -" ,( '™ i) i) (C12) 

= fio ■ zl 1/3 (kaz 3 ) + hi(z) ■ zl 1/3 (kaz 3 ) + /i K (z) ■ zK 1/3 (kaz 3 ). (C13) 
Here we have introduced the following notations (1C81 IC9I) : 
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Ho = J ' dw J r (w)(^2mak 2 w(^ \w)-7i^(w))-k 2 a 2 w 3 ^ 1 \w)^, (C14) 

^(z) = j* dw Tziw) (^2mak 2 w (^ 0) (w) - n^(w)) - k 2 a 2 w 3 ix w (w)\ , (C15) 

Hk(z) = J* dw T K (w)^2mak 2 w (<p (0 \w) - vr (0) H) - k 2 a 2 w 3 n (1) (w)^j . (C16) 
In terms of these definitions we can express the solution </> (1C5I ICIIl ICI2I IC13I) : 



<Kz) = 1 - 21^(1/3) (*f) l ^ zh /3 (kaz 3 ) - 21^(1/3) (^) ^ zK 1/3 (kaz 3 ) 

+ fi zl 1/3 (kaz 3 ) + ni{z)zhp{k(r£) + ^ K {z)zK l/3 {kaz 3 ) + 0{m 2 ). (C17) 

APPENDIX D: THE CANONICAL NORMALIZATION OF THE FIELD AND 
THE PARAMETER L 4 OF THE CHIRAL LAGRANGIAN 

In this section we will present the explicit expressions for the quantities N n and g m ^ intro- 
duced in the section IIIII eq. ( 1331) . 

Let us consider the quadratic part of the 5D action ([6]) where we only take into account the 
first Kaluza-Klein mode of the field <p which corresponds to the pion. We obtain: 

S 5D _> J ' dz (^A 2 R 3 ^lf 2 (z) + ^\d K f}(z)j x jd 4 x \ d,r{x)d^ a {x). (Dl) 



Following ( 1331) let us denote 



. We immediately obtain: 

= (D3) 

where f 2 (z)\ m _ Q may be extracted from (1C17|E2I) . 
The term of the Chiral Lagrangian 

P 4 = tr(D^Z7 + D'*C0tr(C^x + X f ^) = 8-^m^vr a (x)^vr a (x) + C(vr 4 ) (D4) 
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from the AdS/QCD point of view is generated by the action fIDlj) : 



S 5D _, m I d A x d^ a {x)d^ a {x) x / dz ( jU 2 i2V* /J(z)| 



m=0 



+ A 2 i2W ^ + 4" ^/ (^)9 m ^^(z)| m=o ) . (D5) 

If we denote according to (IBHj) 

+ — dzU(z)d m d z f<p(z)\ m=0 I = # m ^2, (D6) 

we obtain: 

U = ^g m ^N~\ (D7) 
where /|(^)| m=0 and S l m /J(js)| m=0 may be extracted from flC7tiCT7l lE2l ) . 

APPENDIX E: THE PARAMETERS L U L 2 AND L 3 OF THE CHIRAL LA- 
GRANGIAN 

In this section we will present the explicit expression for the quantity introduced in the 
section HTT1 eq. (J22D- 

The 4 interaction (the Skyrme-like term in (133]) ) is generated by the quartic part of the 
gauge sector of the full 5D action (1291) : 

S.4 = [ d b x^f abe f cde d tl (t) a d u (t) b d' 1 (t) c d v (t) d . (El) 
J ^9E Z 



The solution ( 1C17I) for <p(z) from subsection ICl is the function f<f,(z) with which we were 
dealing in section IIIIl eq. ( 1301) . We will restrict our consideration to the chiral limit. This 
corresponds to the zero-order solution (f>^\z) (lC5j) : 

- 2^(1/3) J zK 1/3 (kaz 3 )). (E2) 
The normalization factor N<j, is determined, as in (1321) . by the condition (13TT) : 

Zm. A 7 

- m = i. (E3) 
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The effective 4D coupling is (jEIjBDj): 

1/3 



We obtain an effective 4D Lagrangian for the pions: 

^ = g^f^f^d^d^dVd^, (E5) 

7T 

where N n is the normalization factor that determines the proportionality between the fields 7r a 
and (f) a , see the subsection |D| eq. (1D2j) . 
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